Quantum discord for two-qubit CS states: Analytical solution 
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Abstract 

We found an universal local orthogonal transformation which transforms any centrosymmetric 
(CS) density matrix pes into the X density matrix px- H ® HpcsH H = px, where H is the 
Hadamard transformation. Since quantum discord is invariant under the local unitary transfor- 
mations, this remarkable property allows to get the discord of general two-qubit CS states in an 
analytical form using the corresponding well-known formulas for the X states. Examples of systems 
with the CS density matrices are given, including XXZ spin model with the Dzyaloshinsky-Moriya 
interaction, a gas of spin-carrying particles in closed nanopore, and a family of pseudopure states. 
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Introduction. Quantum discord, Q, is a more general measure of quantum correlation than 
entanglement. Discord is important in quantum information theory, quantum metrology, 
condensed matter physics, and so on The evaluation of discord is a very hard problem. 
Even for two-qubit systems, the analytical calculations are restricted to the so-called X 
states. The term "X states" has been introduced in 2007 |2| and denotes the 4x4 density 
matrices which may have non-zero entries only along the main diagonal and anti-diagonal. 
Algebraic properties of such matrices were studied by Rau jsj] . Analytical solution for the 
quantum discord of X density matrices has been obtained in Ref. J4]. On the other hand, 
there is a number of physical systems for which the density matrices have different forms, 
e. g., centrosymmetric ones. The nxn CS matrix is defined by the relations for its elements: 
Qjj = a n-i+i,n-j+i (about the CS matrices see, for example, the reviews [5j and references 
therein). In this communication we establish a connection between CS and X matrices 
through the local orthogonal transformation. Thanks to the fact [l|] that the discord (and 
the entanglement) does not change its value under such transformations, we as a result 
obtain a possibility for the analytic calculation of discord in arbitrary two-qubit CS states. 

CS and X matrices. In quantum mechanics the density matrix must be Hermitian, non- 
negativity defined, and have unit trace. A general 4x4 CS density matrix can be written 
as 
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is the Hadamard transform. The matrix R is orthogonal and R T = R (the superscript 
T denotes a transposition). Taking the matrices (jTJ) ([3j) and performing straightforward 
calculations we establish the following relations 

H ® Hp cs H ®H = Px (5) 

and 

H <g> H Px H ® H = pes- (6) 
Here the parameters of both matrices are related as 
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It is easy to proof the same using the Bloch forms for the density matrices pes and px- 
Indeed, expanding the density matrix ([1]) on the Pauli matrices one obtains 
1, 



Pes = t [1 + ip 4 a x <g> 1 + Ap 2 l <g> a x + 2(p 6 + p 7 )a x ® cr x + 2(p 7 - p 6 )a y <g> cr y 
+ (4pi - l)cr a g) cr 2 - 4p 3 a z ® a y - 4p 5 a y <g> er z ]. 



(9) 



Performing the Hadamard transformations and taking into account that Ha x H = a z , 
HtjyH = —cr y , and Ha z H = a x , we obtain 

H <g> HpcsH ®H = -[1 + 4p 4 0- 2 ® 1 + 4p 2 l <8> o" z + 2(p 6 + p 7 )cr 2 g> a z + 2(p 7 - p 6 )cr y <g> <t b 
+ (4pi - l)a x <&a x + 4p 3 a x ® a y + Ap b a y <g> o x \ = p X - (10) 

The last equality follows from the Bloch form for the matrix (J2]) 

Px = ^{1 - [1 - 2(q 1 + q 2 )}a z <8> 1 - [1 - 2(g x + ?3 )]1 <g> a, + 2(g 4 + g 6 K ® °x 
+ 2(g 6 - g 4 )er„ <g) o y + [1 - 2(g 2 + (lz)\o z ®o z - 2(q 5 - q 7 )a x <g> cr y 
- 2(g 5 + q 7 )a y ® a x } (11) 

and the relations ©. 

As a result, we conclude that the quantum discord of the state pes is expressed through 
the discord of the state px- 

Q(pcs) = Q(px), (12) 

where the entries of p x are given by Eqs. ([7j). Together with the analytical formulas for the 
discord of X states |4j , this completes our solution. 

Physical examples. As a first illustration, consider the anisotropic XXZ model with the 
Dzyaloshinsky-Moriya interaction. When the Dzyaloshinsky vector D is oriented along re- 
direction, the Hamiltonian for the two-qubit chain reads [6] 



U = Jo*o* + Ja\al + J z a{a z 2 + D x {a\a z 2 - o\o\). 



(13) 



Here J and J z are the coupling constants and af (i — 1,2 and a = x,y,z) are the Pauli 
matrices. In open form the Hamiltonian is given as 
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This matrix is centrosymmetric. Because the sums and products of CS matrices are again 
the CS matrix, the corresponding Gibbs density matrix is CS one. Therefore, the thermal 
quantum discord is found here in an exact analytical form. 

Another example is related to the dynamics in NMR of quantum correlation (discord) for 
the pair of nuclear spins in a nanopore filled with a gas of spin-carrying molecules or atoms 
Q]. The corresponding reduced density matrix is also the CS one: 
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In these relations, iV is the number of particles confined in a nanopore, a is the normalized 
coupling constant, and (3 is the inverse dimensionless temperature. In the paper it was 
succeeded to calculate the quantum discord only for a particular case p = u = and q = r 
when the density matrix (fl~5l) is reduced to the Bell-diagonal form. Applying the method 
developed above we see that after performing the Hadamard transformation 03), the matrix 
(1T5]) takes the X structure 
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or in the Bloch form: 



[1 + 2p(a z 1 + a z 2 ) + Ara\ai + Aqa\a z 2 - u(a^ + a\a x 2 )\. 



(18) 
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Using the general formulas of Ref. 4 we can in principle get the discord. But there is a 
simpler way. Indeed, perform an additional local unitary transformation to eliminate 'xy' 
cross-terms in Eq. f ll8p , that is to reduce the density matrix (1T7|) to the real X form. We are 
achieving this goal with the transformation U ®U p'U + ®U + = p", where U = exp(—i<pcr z /2) 
and (f = — | arctan(2w/r). After this transformation, the density matrix (j!7p takes the form 
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Jsing now the formulas for calculating the quantum discord for the real X density matrices 
8] we finally obtain 

g = min{Q 1 ,Q 2 }, (21) 
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In Eqs. fl22|) and ( 1231 . 5 and SV are the entropies of the full and reduced density matrices, 
respectively: 

S = -]TA,log 2 A„ (25) 
where the eigenvalues Aj of the density matrix under consideration are given as 



Ai >2 = - + q± [p 2 + (2wsin2y? - r cos 2^) 2 ] 1/2 , 



A 3 ,4 = --q±\r\ 
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FIG. 1: Quantum discord in nanopore as a function of dimensionless time. 

and 

Sr = ~{\ + P) log 2 (^ +p)-{\-p) log 2 (i - p). (27) 
Typical time dependence for the pairwise quantum correlation of spin-carrying particles in 
a nanopore is shown in Fig. [TJ The correlation oscillates between zero and the saturation 
value which is Q when p = u = and r = q = |tanh 2 (/V2) [see Eq. (26) in Ref. 7]. 
Our third example concerns with the pseudopure (PP) states 

1 - a 

ppp = a\^)^\ + —I J (28) 

where is an arbitrary two-qubit pure state, I is the identity operator, and the proba- 
bility a G [0, 1]. The states like (I2"8"j) are studied as a possible resource for NMR quantum 
computing (see, e. g., and references therein). It is easy to check that if 

|^> = o(|00> + |ll»+6(|01> + |10» (29) 

(|a| 2 + \b\ 2 = 1/2 and {|00), |01), |10), |11)} is the computational basis) the state ppp will be 
CS one and we can calculate its quantum discord using the above scheme. 

Conclusions. We have established the relation between the CS and X matrices via the 
universal local orthogonal transformation. This allows to find the discord of any CS states 
using available formulas for the discord of X states. CS quantum states appear in different 
important physical problems, three of which have been discussed. 
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